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ABSTRACT 

We use signal enhancement techniques and a matched hlter analysis to search 
for the K band spectroscopic absorption signature of the close orbiting extrasolar 
giant planet, HD 189733b. With timeseries observations taken with NIRSPEC at the 
Keck II telescope, we investigate the relative abundances of H2O and carbon bearing 
molecules, which have now been identified in the dayside spectrum of HD 189733b. 
We detect a candidate planet signature with a low level of significance, close to the 
~ 153 kms^ 1 velocity amplitude of HD 189733b. However, some systematic variations, 
mainly due to imperfect telluric line removal, remain in the residual spectral timeseries 
in which we search for the planetary signal. Using principal components analysis, the 
effects of this pattern noise may be reduced. Since a balance between the optimum 
systematic noise removal and minimum planetary signal attenuation must be struck, 
we find that residuals, which are able to give rise to candidate planet signatures, 
remain. The robustness of our candidate signature is therefore assessed, enabling us 
to conclude that it is not possible to confirm the presence of any planetary signal 
which appears at F p /F* contrasts deeper than the 95.4 per cent confidence level. Our 
search does not enable us to detect the planet at a contrast ratio of F p / F* = 1/1920 
with 99.9 per cent confidence. 

Finally, we investigate the effect of model uncertainties on our ability to reliably 
recover a planetary signal. The use of incorrect temperature, model opacity wave- 
lengths and model temperature-pressure profiles have important consequences for the 
least squares deconvolution procedure that we use to boost the S/N ratio in our 
spectral timeseries observations. We find that mismatches between the empirical and 
model planetary spectrum may weaken the significance of a detection by ~ 30 - 60 per 
cent, thereby potentially impairing our ability to recover a planetary signal with high 
confidence. 

Key words: Line: profiles - Methods: data analysis - Techniques: spectroscopic - 
Stars: late-type - Stars: individual: HD 189733 - Stars: planetary systems 



1 INTRODUCTION 

The field of exoplanet spectroscopy has advanced rapidly 
since the first observations of secondary eclipse events en- 
abled the bri ghtness temperatures of transiting planets to b e 
determined l|Charbonneau et ai]|2005l : IDeming et al.ll2005l) . 

* E-mail: j.r.barnes@herts.ac.uk 



thereby confirming heating due to stellar irradiation. Thanks 
largely to the Spitzer Space Teles cope (|Houck et al.l |2003 : 
iFazio et alj|2004l ; iRieke et alj|2004 ). the spectral energy dis- 
trib ution of a num b er sys tems has now been estimated (e.g. 
see iBurrows et alj (|2008h ). It has recently been suggested 
that close orbiting extrasolar giant planets (CEGPs) can 
be divided into two possible sub-groups, one which exhibits 
absorption spectra and one which exhibits emission features 
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due to the presence of a stratosphere l|Burrows et alj [20081 ; 
iFortnev et al.l [2008) . Individ ual planetary atmosph eres are 
no doubt more complex with iBurrows et all l|200Sft . for ex- 
ample, preferring to use parameterisations of the degree of 
stratospheric absorption and heat redistribution. Neverthe- 
less, H P 209458b appears to fall broadly into the latter cat- 
egory (jKnutson et al.l I2OO8I ) while the spectral energy dis- 
tribution of HD 189733b is f ound to be we ll fit by models 
where no stratosphere forms (|BarmaJ2QOSh . 

Although planetary signatures are easier to detect at 
mid-infrared wavelengths, characterisation of atmospheres 
is easier at shorter wavelengths where observational sensi- 
tivities are sufficient to probe the large amplitude signatures 
of species such as H2O and carbon bearing molecules. The 
first evidence for water and organic molecules in the atmo- 
spheres of CEGPs has come from space based broadband 
photometric obs ervations and l ow res olution transi t spec - 
troscopy. While iBarmanl (|2007h and iTinetti et ail (|2007l ) 
claimed independent dete ctions of H2O at different wave- 
lengths, [SwahTetlil] (|2008l ) have now identified CO and CH4 
in transit spe ctra which probe the terminator of the planet 
HD 189733b. iGrillmair etaD (|2008l ) reported a detection of 
water in the dayside spectrum through Spitzer space tele- 
scope observations of secondary eclipse events below 7.5 \im. 
Most recently NICMOS/HST data have been used to infer 
the additional prese nce of CO and CO 2 in the dayside spec- 
tra of HD 189733b |Swain et al.ll2009l . hereafter S09). 

Rather than attempting to identify molecules from their 
broadband spectral signatures, we present a method which 
attempts to identify spectral structure through a statisti- 
cal examination of the many thousand individual transitions 
found in a typical high resolution (R ~ 25,000 - 50,000) plan- 
etary spectrum. Specifically, in order to detect the faint plan- 
etary signal, spectral deconvolution is applied to individual 
spectra in a timeseries in order to derive mean absorption 
profiles with boosted S/N ratios. Modelling the phase de- 
pendent contrast ratio and radial velocity motion enables 
the maximum planet/star contrast ratio and velocity am- 
plitude to be measured (from which the true mass of the 
planet may also be determined). Crucially, since the planet 
need not be transiting, this spectroscopic method increases 
the sample of planets which can potentially be studied with 
current instrumentation. With the aim of extending space 
based mid-infrared measurements of planet / star contrast ra- 
tios into the near infrared, we search for the CEGP signa- 
ture, manifested as H2O, CO and CO2 absorption in high 
resolution K band spectroscopic timeseries observations of 
HD 189733 (K1V-K2V). 



2 DATA REDUCTION 
2.1 Observations 

K ban d observations of HD 189733 were secured with NIR- 
SPEC (|McLean et al.lll998l ) at the Keck II Telescope on UT 
2008 June 15 and June 22. Respectively, a total of 219 and 
154 spectra were recorded using a 1024 x 1024 InSb Aladdin- 

3 array. With the NIRSPEC-7 blocking filter, a wavelength 
span of 2.0311 ^.m- 2.3809 urn was achieved with a slit width 
of 0.432", giving a resolution of R ~ 25, 000. Our 60 sec 
exposures comprised of 12 coadds, each of 5 sees duration. 




Figure 1. HD 189733b orbital phase diagram. Observations made 
on June 15 cover phases tp = 0.303 — 0.429 while observations 
made on June 22 cover <j> = 0.498 — 0.581. Observations in the 
range <f> = 0.498 — 0.517 were not used in the subsequent analysis 
since the planet is eclipsed during these phases. Only the phases 
used for analysis are shown. 

The observations are summarised in Table [1] The seeing was 
mostly good at around 0.6-0.7" although observations were 
plagued by cloud for a period on June 15. 

2.2 Data extraction 

Pixel to pixel variations were corrected for each frame using 
flat-field exposures taken with an internal tungsten refer- 
ence lamp. The worst cosmic ray events were removed at 
the pre-e xtraction stage u sing the Starlink FIGARO routine 
BCLEAN (|Short r idge 1993). Since we chose not to use an 
ABBA nodding sequence, in order to maximise stability 
in the K band, we carried out the same extraction proce- 
dure as detailed in iBarnes et ail (l2007bl ) for previous ob- 
servations of HD 189733. The spectra were extracted using 
echomop's implem entation of th e optimal extraction algo- 
rithm developed bv lHorne] (1 19861). echomo p rejects all but 
the strongest sky lines IBarnes et al.l (|2007bl ) and propagates 
error information based on photon statistics and readout 
noise t hroughout the extract ion process. 

In IBarnes et al.l (|2007bl ). we reported observations of 
HD 189733 which were made in conditions of variable cloud 
and seeing. This resulted in a spectral timeseries in which 
the S/N ratio of individual spectra varied greatly. We found 
that rejecting the lowest S/N frames increased the sensitiv- 
ity of our method. We thus carried out a similar rejection 
procedure, set at an arbitrary S/N = 350 to reject the low- 
est S/N ratio frames which formed a distinct distribution 
separate from those frames which were observed in the best 
seeing conditions. A total of 63 frames were rejected from 
observations made on June 15 (S/N ratios of the rejected 
frames were 245 ± 127) while none were rejected from ob- 
servations made on June 22; a clear reflection of the more 
variable conditions prevailing on the first night. The S/N 
ratios of the residual timeseries were 497 ±125 and 570 ± 77 
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Table 1. Keck/NIRSPEC observations of HD 189733 for UT 2008 June 15 and 22. 



3 



UT Date 


UT start of 


UT start of 


Time per 


Number of 


Number of 




first frame 


last frame 


exposure [sees] 


co-adds per frame 


observations 


2008 June 15 


08:36:29 


15:18:33 


5 


12 


219 


2008 June 22 


10:43:37 


15:08:46 


5 


12 


154 



on June 15 and 22 respectively. The phases of observations 
used during the subsequent analysis stages are represented 
visually in Fig. 1. 



2.3 Residual spectra 

Follo wing the pro c edure s outlined in iBarnes et al l l|2007bl ) 
and IBarnes et al.l ((2008) , we attempt to extract the plan- 
etary signature from timeseries spectra by removing the 
dominant spectral contributions; namely the stellar spec- 
trum and the telluric lines. A master spectrum is gener- 
ated by combining all observed frames on a given night af- 
ter nearest-pixel alignment of each individual order in each 
spectrum to minimise blurring. The master template is sub- 
tracted from each spectrum in turn after shifting, scaling 
and blurring/sharpening (again on an order to order basis). 
This latter procedure is achieved by calculating the first to 
fourth order derivatives of the template spectrum. The mas- 
ter template is then scaled to each observed spectrum by us- 
ing a spline with a fixed number of knots. Scale factors for 
the derivatives are similarly calculated. A model of each ob- 
served spectrum can thus be calculated by means of a Taylor 
expansion of the template spectrum . The procedure is de- 
scribe d in detail in Appendix A of ICollier Cameron et al.l 
(2002). Since we observe HD 189733b at phases close to 
phase <f> = 0.5 when the planet spectrum is Doppler shifted 
on the steepest part of the radial velocity curve, the mas- 
ter template spectrum contains only a very weak, blurred 
out copy of the planetary spectrum. Hence subtraction of 
the template does not significantly attenuate the planetary 
signature. 



2.4 Residual pattern noise removal 

The resulting residual timeseries should thus ideally only 
contain noise and a copy of the planetary spectrum. How- 
ever, residual pattern noise remains in the subtracted time- 
series. Examination indicates that the dominant patterns 
are time variant residuals which are introduced when sub- 
tracting the scaled master template during the procedure 
outlined in §2.3. Although this procedure is effective at re- 
moving the stellar and telluric lines, variations between the 
strength of these lines throughout the night (due to chang- 
ing airmass and hence telluric strength) are sufficient to pre- 
clude consistent results over timescales of the variations. 

In addition to pattern noise induced by telluric ef- 
fects, no ise may arise from t ime variable fringing effects, as 
found bv lBrown et~al] (|2002h for K band NIRS PEC observa- 
tions for example. To remove these two effects, iBrown et al.l 
l|2002i ) carried out a two stage procedure in the form of a 



regression and singular value decomposition filtering pro- 
cedure on their spectra. Our implementation of the Tay- 
lor expansion algorithm, me ntioned above and in detail in 
ICollier Cameron et al.l ((2002), is analogous to the regression 
step. In a similar proce dure to the singular value decompo- 
sition implemented by IBrown et al.l i|2002l ), the remaining 
pattern nose can be removed by calculating the correlation 
matrix of the time variations in individual spectral bins for 
our set of residual spectra. The eigenvectors of the corre- 
lation matrix, which account for the largest fraction of the 
variance, are the principal components of the residual spec- 
tra. By subtracting only the principal components, any re- 
mai ning pattern noise can be furt her removed (see Appendix 
B of lCollier Cameron et al.ll2002l for further details). A bal- 
ance must be struck between removing pattern noise and 
maintaining information in the residual spectra. We found 
that shifts of 1-2 pixels during a typical night of observations 
could be attributed to points during the night where the 
instrumental configuration was slightly modified to enable 
observations of a star for another project. We thus applied 
principal components analysis to continuous blocks (i.e. be- 
tween observations of the other star, where the 1-2 pixel 
shifts occurred) of HD 189733 observations and found that 
2 to 3 components were necessary to remove the remaining 
pattern noise in the data without attenuating the planetary 
signal (see §4). 



3 DECONVOLUTION AND PLANETARY 
MODELS 

In order to extract the planetary signature from the cor- 
rected residual spec tra, we used least squares deconvolution 
IIDonati et al. 1997) which requires the use of a model spec- 
trum to describe the strengths (normalised profile depths) 
and wavelength positions of the strongest planetary opaci- 
ties over the wavelength range o f our observations. Our im- 
plementation of the algorithm (|Barnes et al.l [l998h propa- 
gates errors from the input spectra a nd has been used in re- 
flected light searc hes in the opt i cal by Collier Cameron et al.l 
i| 19991 . 12OO2I I and iLeigh et all (|2003al lbh. For each residual 
spectrum a deconvolved profile is obtained, potentially con- 
taining a copy of the Doppler shifted planetary profile which 
can then be detected owing to the effective boost in S/N ra- 
tio. Typical boosts in S /N ratio of a few times to a few tens 
of times are achieved since several hundred to several thou- 
sand planetary absorption lines are used to deconvolve the 
planetary signature. The 2.04- 2.06 \im and 2.36- 2.38 wave- 
length ranges are omitted in all the following analyses due 
to the strong telluric features which dominate these regions 
of the spectra. 

We have generated several models to represent 
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the emergent spectrum of HD 189733b. Our standard 
HD 189733b m odel, an updated ve rsion of the model de- 
tailed in §3.2 of[B arncs et all (|2007bh . is generated for an at- 
mosphere with solar metallicity, a temperature, T = 1250 K 
and surface gravity, log g = 1.33 ms _1 . For a detailed de- 



script io n of the model o p acities and se tup seelFerguson et al.l 
l|2005h . iBarman et all <|200ll ) and iBarman et al.l (120051 1 
BHA05. The most recent models and their ability to fit 
HD 189733b observations made with the Spitzer Space Tele- 
scope are discussed at length in IBarman! (|2008h . A number 
of further models with adjusted temperature pressure pro- 
files and relative chemical abundances will be discussed in 
the following sections. The S/N ratio after deconvolution 
with the standard model yields timeseries with mean profile 
S/N ratios of 7520 and 9450 on June 15 and 22 respectively. 
The mean deconvolved wavelength of the spectra depends 
on the wavelength dependent count rate and the relative 
strengths of the absorption lines. For the standard model, 
Xmean = 2.19 urn. As reported in §4, the mean wavelength 
may change slightly when the line list used for deconvolution 
is modified. 



4 RESULTS 

A matched Gaussian filter approach is used to search for 
the planetary signal, where a model describing the radial 
velocity shift and planet/star contrast phase function is 
used to search for the deconvolve d absorption profile of 
the planet in the spectral ti meseri es (|Collier Cameron et al.l 
2002|; iBarnes et aLI l2007al lbl. [ 20081 ). Pairs of maximum con- 
trast ratio (eo) vs velocity amplitude (K p ) are used in a two 
dimensional \ 2 search to find the combination which yields 
the best improvement in x 2 ■ The significance of candidate 
enhancements in x 2 are assessed by a bootstrap procedure 
which randomis es the order of the data wit hin each night 
of observations l|Collier Cameron et alj|2002h . This process, 
carried out several thousand times, scrambles any planetary 
signature, but enables the data to retain the ability to give 
false \ 2 enhancements due to systematics which may remain 
in the data above the photon noise. Reliable confidence lev- 
els can thus be plotted on the 2-parameter x 2 landscapes of 
logio(eo) vs (K p ) x 2 - 

Calibration of contrast ratios is achieved by injecting a 
"fake" model planet spectrum into the timeseries and then 
recovering the signal using our matched filter method. The 
fake planetary spectrum is injected with known eo and K p 
after extraction of the spectra and before any of the subse- 
quent steps described above are carried out. In this way, we 
are also able to a ssess our ability to correct ly recover a plane- 
tary signature (|Barnes et al.ll2007bl . [2008t ). We find that for 
a planet recovered with high significance, there may be a 
slight shift in the recovered K p velocity. This is most likely 
due to some removal of planetary signature when subtract- 
ing the template star and during principal components anal- 
ysis. During these procedures, we chose parameters which 
strike a balance between removing residuals while not signif- 
icantly affecting the planet signal. Essentially, an absorption 
signature located at phases close to tj> = 0.25 will be attenu- 
ated most, since, in this region the planet shows the small- 
est radial velocity gradient with orbital phase. The resulting 
magnitude of the velocity amplitude uncertainty is typically 



< 5 km s _1 for a planet simulated with greater than 99.9 
per cent significance but may be as much as 10- 20 km s _1 
where a fake planetary signature is injected with ~ 95.4 per 
cent significance. This is especially true if the data contain 
systematics above the photon noise level of the data. 

4.1 Standard model 

In Fig. [2] (top left) we present the phased deconvolved time- 
series of the residual spectra (i.e. spectra with removed stel- 
lar spectrum and tellurics and containing a potential plan- 
etary signature) based on our standard model (see §3). For 
plotting purposes only, the timeseries has been normalised 
using the formal variances since some phases (particularly 
4> ~ 0.372 — 0.429 at the end of the first night of observa- 
tions ) are more noisy than others. This enables the noise 
structure at all phases to be more clearly seen. The black 
and white greyscale values are set at ±1.5<r in the plotted 
normalised timeseries. It should be stressed that the true for- 
mal variances are utilised when searching for the planetary 
signal in the un-normalised deconvolved timeseries. Hence 
spectra with lower S/N receive a lower weighting in our 
analysis. Since the planet undergoes eclipse during phases 
(f> = 0.498 — 0.517, we do no use spectra taken during this 
interval in our analysis. These spectra are however plotted 
for completeness in Fig. &[2] 

In Fig. [2] (top right) we present the x 2 landscape plot of 
log(eo) vs K p based on our standard model. Dark features in 
the plot represent enhancements in x 2 whose significance can 
be measured relative to the plotted confidence levels. The 
large black feature representing the greatest enhancement in 
X 2 appears with low confidence (in the 68.3 - 95.4 per cent 
confidence region) at log(eo) = -3.41 and K p = 85 kms . 
However, as can be seen in the phased timeseries, a number 
of low-level residual features are present. These appear as 
dark absorption areas, covering localised regions of velocity 
and phase. We believe that these features are responsible for 
the K p — 85 kms -1 signature and result from imperfect re- 
moval of telluric and stellar lines during our analysis, giving 
rise to false signals . As has been demonstrated previously 
l|Barnes et al.ll20 08h. a clear detection of the planetary sig- 
nature would be expected to result in a more localised x 2 en- 
hancement and greater significance than the K p = 85 kms -1 
X 2 enhancement. Nevertheless, K p is known for HD 189733b 
(since the system is eclipsing and the orbital inclination is 
known) and is indicated by the dashed lines in Fig. [21 A can- 
didate planetary signature should thus appear at, or close 
to (see above), this velocity amplitude. 

For HD 189733b, the standard model predicts a maxi- 
mum contrast ratio of logio(eo) = -3.163 or F p /F* — 1/1460 
over the wavelength span of our observations. We are how- 
ever unable to detect the planetary signature, at the mean 
deconvolved wavelength of 2.19 um, with 68.3, 95.4, 99 
and 99.9 percent confidence levels of logio(eo) = -4.065, - 
3.491, -3.366 & -3.193 or F p /F, = 1/11600, 1/3100, 1/2320 
& 1/1560 respectively. In light of much better observing con- 
ditions, this is a si gnificantly more sens itive result than our 
2006 observations IjBarnes et alj|2007bh permitted. Consid- 
erable care must be exercised if quoting sensitivities at con- 
trasts ratios deeper than the 95.4 per cent level (this is in- 
vestigated further in section §4.3 below) owing to candidate 
signatures which arise from systematics at these levels. HD 
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Figure 2. Phased deconvolved timeseries spectra of HD 189733b and corresponding 2-dimensional \ 2 landscape plots for analysis using 
the standard model (top) and enhanced CO2 model (bottom). Left: The dashed line in the timeseries plots represents the position of 
the planet as a function of phase. The black and white levels are set at the 1.5ct level of the timeseries. Right: \ 2 plot for matched filter 
combinations of maximum contrast ratio, logio(eo) vs K p . Observations covering phases <f> = 0.498 — 0.517 are shown for completeness 
but were not used in the analysis since the planet undergoes eclipse at these phases. Black and white represent the best and worst 



spectra. No enhancement in \ 2 is measured at the known velocity amplitude, K p = 152.6 km s _1 , marked by the dashed line. Upper 
confidence levels of 63.8, 95.4, 99 &l 99.9 per cent (top to bottom solid lines) are plotted. The bold + symbols mark the mean value of 
the contrast ratio over the range of observations for the standard model (top right) and the corresponding mean value of the contrast 
ratio from the NICMOS/HST S09 observations (bottom) and arc logio(eo) = -3.163 & -3.286 (F p /F* = 1/1460 & 1/1930) respectively. 
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Figure 3. Comparison of the standard model (black) and 
augmented CO2 abundance model (grey) with the observed 
NICOMS/HST dayside spectrum of S09. The enhanced CO2 
model shows a considerably improved fit in the 1.9- 2.2 um region 
relative to the standard model. 

189733b is therefore not detected at K p = 152.6 km s _1 at a 
contrast which is 2.1 and 1.1 times deeper (95.4 per and 99.9 
per cent confidence respectively) than the standard model 
predicts. 

4.2 Enhanced CO2 model 

Recent Spitzer/IRAC l|Grillmair et all I2008T ) and 
HST/NICMOS (S09) observations have enabled the 
dayside spectrum of HD 189733b to be measured using low 
resolution spectroscopy. In order to reliably fit the spectrum, 
a greater than expected abundance of CO2 is required with 
S09 reporting CO2 mixing ratios of c ~ 0.1 — 1 x 1CT 6 . We 
have generated a model with augmented CO2 abundance 
which enables us to match the HST/NICMOS observations 
(Fig. 3). From here on, in §4, we only consider this model. 
We emphasise that our current model does not contain 
some of the hot CO2 bands which have been identified 
by Fourier Transform Spectroscopy carried out at the Jet 
Propulsion Laborat ory and include d in t he latest edition 
of HITRAN (see iRothman etafl 120091 and references 
therein). We therefore caution that in order to achieve the 
required level of absorption from CO2 in the 1.9-2.2 urn, 
the relative strengths of individual opacities in the output 
model are likely overestimated. Nevertheless, based on the 
S09 results, inclusion of such opacities would appear to give 
a more accurate representation of the expected opacities 
found in the spectrum of HD 189733b. Fig. 2 (bottom 
left and bottom right) presents the deconvolved timeseries 
and log(eo) vs K p \ 2 plot after deconvolution using our 
augmented CO2 model spectrum. At K p — 152.6 kms -1 , 
the expected planetary signature is not detected with 
68.3, 95.4, 99 and 99.9 per cent confidence levels of 
login (e ) = -4.074, -3.529 & -3.400 & -3.283 or F p /F* = 
1/11600, 1/3380, 1/2510 & 1/1920 respectively. The sen- 
sitivities are slightly greater than for the standard model, 
although we note that the enhanced CO2 sensitivities are 
quoted for a centroidal wavelength of 2.15 urn rather than 
2.19 urn. This shift in mean wavelength results from greater 
normalised depths of the enhanced CO2 at the shorter 
wavelengths of our observations. In addition, although 



the mean planetary flux level is lower in the regions with 
enhanced CO2, the recorded count rate in the observed 
star + planet spectra is higher, leading to higher contrast 
confidence limits with this model. The mean 2.0-2.4 urn 
planet/star flux ratio reported by S09 is logio(eo) = -3.286 
or F p /F, — 1/1930 (i.e. almost identical to our 99.9 per 
cent confidence level of 1/1920), indicating that we are 
sensitive down to the 99.9 per cent level. 

4.3 Wavelength splitting the data - a candidate 
signature? 

While there is no clear candidate signature at the expected 
velocity amplitude of the planet in the enhanced CO2 de- 
convolved timeseries, we do detect a signal with relatively 
low confidence (95.4 percent) at K p — 167.9 km s -1 and 
Iogio(eo) = -3.63 [F v /F* = 1/4270). Although it is likely 
that this signature is again the result of alignment of sys- 
tematic absorption features in the timeseries at this ve- 
locity amplitude, we have investigated splitting the time- 
series data into two spectral regions. Deconvolution was car- 
ried out on the first three orders (region 1: 2.03-2.18 u.m) 
and on the second three orders (region 2: 2.21 - 2.36 u.m) 
independently before carrying out a search for the plan- 
etary signature. With the data split in this manner, re- 
gion 1 contains H2O and CO2 opacities while region 2 
contains H2O and CO (bandhead at ~ 2.29 urn) opac- 
ities. While region 1 did not reveal any candidate sig- 
nature close to K p — 152.6 kms -1 , region 2 has a more 
well defined candidate signature at K p — 147.8 kms -1 with 
logio(eo) =-3.449 (1/2810) and 97.2 per cent confidence. 
The S09 contrast ratio for HD 189733b over the wavelength 
span of region 2 is (F p /F* ~ 1/1340). To assess the true 
nature of the region 2 signature, we refer the reader to 
Fig. 4 (left) which again indicates that there are a number of 
systematic features. Close examination reveals that the ex- 
pected velocity position of the planet as a function of phase 
(indicated by the dashed line) appears to pass through, or 
near, a number of contiguous absorption regions. To inves- 
tigate the contribution of these regions to the candidate sig- 
nature, we have carried out three tests. 

1) Analysis of the data observed on 15th June - 
4> — 0.303 — 0.429 alone. The result is a candi- 
date signature with K p = 153.9 kms -1 and with 
logio(eo) = -3.008 (1/1020) and 98.8 per cent confi- 
dence. 

2) Analysis of the data observed on 22nd June - 
4> = 0.517 — 0.581 alone. A candidate signature with 
K p = 141.8 kms" 1 and with logi (e ) = -3.528 (1/3372) 
and 94.8 per cent confidence is found. 

3) Contiguous dark regions omitted from the analysis by 
eye (the regions are not of sufficient amplitude to enable 
reliable sigma-clipping). No candidate signature within 
26 kms -1 of the known K p — 152.6 kms -1 is apparent. 

The varying velocity and contrast ratio of the candi- 
date signals from tests 1 and 2 suggest that if any plan- 
etary signature contributes to the \ 2 enhancements, it is 
biased by some other factor. Test 3, in which the contigu- 
ous regions are omitted from the analysis, has the effect 
of removing the candidate signature seen in Fig. 4 (right) 
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Figure 4. As for Fig. 2, but using only the wavelength range 2.21-2.36 |om. The corresponding NICMOS/HST S09 contrast ratio for 
this region is marked by a + symbol with F p /F* ~ 1/1340. A candidate signature is detected with 97.2 per cent confidence close to 
the expected K p velocity, although the contrast ratio of logio(eo) = -3.449 (F p /F t = 1/2810) is lower than expected. See main text for 
details. 



completely. The confidence levels, at K p — 152.6 km s , 
with the omitted contiguous residual absorption regions 
are log 10 (e ) = -3.964, -3.355 & -3.214 & -3.098 or F p /F, = 
1/9200,1/2260,1/1640 & 1/1250 respectively. In other 
words, the phases which do not show contiguous blocks of 
absorption residuals in the timeseries (66 per cent of the 
recorded spectra) along the radial velocity path of the planet 
do not possess the ability to recover a planetary signal. 

We stress that the argument asserting that the contigu- 
ous residual absorption regions are wholly due to systematics 
and solely responsible for producing candidate signatures is 
however not strictly true. Any residual absorption features 
in the timeseries have the ability to modify the contrast ra- 
tio and velocity amplitude of a true planetary signal. Since 
the residual absorption features may be expected to vary 
in strength it is not unlikely that they would result in a 
planetary signature modified by differing degrees in tests 1 
and 2. In our third test, removing 33 per cent of the data 
along the expected radial velocity curve of the planet leaves 
only regions which are consistent with the mean level, or 
regions of contiguous "emission" relative to the mean level. 
One might expect that this procedure would severely im- 
pair our ability to detect a planetary absorption signature. 
In the hypothesis that dark regions are artifacts of the data 
processing (i.e. imperfect telluric/stellar line removal), after 
their removal (test 3) we can rule out our ability to detect 
the planetary signal with 95.4 per cent confidence at a level 
of F p /F* — 1/2260. In light of this and our 99.9 per cent up- 
per limit (§4.2) on the contrast ratio, we believe that since 



S09 detect the planet with F p /F, = 1/1930, further investi- 
gations of model dependency on our analysis are required. 



5 MODEL DEPENDENCY EFFECTS 

We are confident that the planetary signature is not severely 
attenuated (there is inevitably some attenuation as de- 
scribed in §4) during our analysis procedure since fake plane- 
tary signals which are injected before analysis are recovered. 
A cause of our inability to detect a planetary signature is 
likely to stem from a mismatch between the model plane- 
tary spectrum and observed planetary spectrum. The most 
likely direct causes of line strength mismatch a nd model line 
wavelen gth opacities were first highlighted in iBarnes et al.l 
l|2007al ). Line strength mismatches may arise from incorrect 
treatment of the model atmosphere, including uncertainty in 
the exact form of the temperature-pressure (T-P) profile. In 
addition, the precision of the calculated opacities is limited 
by the accuracy of the Einstein A coefficients. This latter 
effect may be true for im portant molecular species such as 
H2O (|Barber et al.l l2006) for instance. We investigate rela- 
tive line depth, temperature and wavelength uncertainties 
below. 



5.1 Relative line depths 

Although the model planetary spectrum may show little 
variation as a result of T-P profile changes when observed 
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at low resolution, the relative line depths may change sig- 
nificantly. In addition, the gradient of the T-P profile de- 
termines the absolute strength of the absorption lines. To 
investigate these effects, we have generated a series of ad 
hoc T-P profiles and resulting emergent spectra. Fig. 5 
shows the models plotted for a short region of wavelength 
space. Steeper T-P gradients lead to the formation of deeper 
lines whereas the relative line strengths vary from model to 
model. These effects are important since mismatch of the 
model and observed spectra line depths will lead to a non- 
optimally deconvolved line profile, and therefore decrease in 
sensitivity. If all the lines are less deep, they will simply be 
harder to detect above the noise level. 

To investigate the degree to which our ad hoc models 
affect the sensitivity of our procedure we used our stan- 
dard model (black in Fig. 5) to inject a fake planet into the 
HD 189733 timeseries. The fake planet was then recovered 
and calibrated to match the contrast ratio at which it was 
injected. By deconvolving with line lists derived from each 
of the different models shown in Fig. 5, we find that we 
are able to recover the planetary signature in all cases. The 
planetary signature is however recovered with an incorrect 
contrast ratio and modified relative confidence. Fig. 6 shows 
the relative contrast ratio for a fake HD 189733b planet in- 
jected into the data with 99.9 per cent significance. Model 
represents the standard model calibration to which the 
simulations are normalised. The contrast ratio is incorrectly 
recovered, with model 1 showing a 2.3 per cent overestima- 
tion of the contrast level and model 4 indicating a 66 per 
cent underestimation of the contrast ratio. In all cases how- 
ever, the ad hoc models appear to recover the planet with 
increased significance. The effect is nevertheless relatively 
small, with models 2 and 3 showing the greatest increase in 
confidence. Model 2 indicates an increase in significance of 
12.5 per cent relative to the 99 per cent confidence level. One 
might naturally expect a decrease in confidence to arise from 
mismatch of the line strengths during deconvolution rather 
than the counter-intuitive increase. We believe that the in- 
crease is most likely due to models 1-4 yielding strong lines 
which become even stronger and weak lines which become 
weaker relative to the standard model. If one of these models 
were a closer match to the empirical HD 189733b spectrum, 
we note that the relative significance of the standard model 
would decrease (with a maximum reduction in sensitivity of 
14.8 per cent relative to the 99 per cent confidence level). In 
conclusion, the above ad hoc models alone are not able to 
explain the lack of true planetary signal (our observations 
are after all still sensitive enough to detect HD 189733b) 
through mismatch of line strengths. Relative to the con- 
fidence levels in Fig. 6, a true planetary signal could not 
change its confidence by more than half of the separation of 
the 99 and 99.9 per cent confidence levels. 

5.2 Opacity wavelength uncertainties 

Moderate wavelength uncertainties lead to an effective 
degradation of deconvolved resolution while model temper- 
ature uncertainties may lead to line strength mismatches 
with the observed spectra during deconvolution. We were 
able to use an unp ublished improved version of the BT2 
l|Barber et al.ll2006l ) water line list (with reduced wavelength 
uncertainties) to investigate these effects. Most of the strong 
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Figure 5. Ad hoc T-P profiles and resulting spectra for 2.135- 
2. 145 um region. The black profile and spectrum represents the 
standard, standard model. While a steeper temperature profile 
at pressures where lines predominantly form (e.g. red/blue mod- 
els) result in spectra with deeper absorption lines, note that the 
relative line strengths also differ from model to model. 



lines which contribute to the deconvolution are transitions 
between states whose energies are experimentally known to 
very high accuracies. Consequently, by selecting only those 
lines which are greater in strength than 1/10,000 of the 
strongest line, we eliminate a large number of lines which 
are expected to have larger uncertainties in calculated po- 
sitions and which in any case have negligible contribution 
to the deconvolved profile. Where the energies of the up- 
per and lower states of a transition are both known experi- 
mentally, these are used for modified BT2 line frequencies, 
rather than the ab initio calculated values. At 1250 K, 75 
per cent of the water lines in our trimmed list are transi- 
tions between experimentally-known levels, and only 25 per 
cent employ BT2 ab initio frequencies/ wavelengths (in all 
cases, however, BT2 Einstein A values are used in comput- 
ing line strengths as these are generally more accurate than 
exp erimentally-determin ed values). 

iBarber et al.l (|2006h state that a comparison of the BT2 
ab initio frequencies with experimentally-known transitions 
shows that the positions of ~ 40 per cent of the lines tested 
are accurate to within 0.1 cm -1 and 91 per cent are within 
0.3 cm -1 . At 2.2 [im, this corresponds to resolutions of R = 
45,500 and R = 15,150 respectively. Clearly, since water is 
the dominating opacity, these uncertainties will play a role 
in degrading the resolution of a deconvolved profile for data 
sets with resolutions of R^ 15,000. These uncertainties are 
therefore applied to the 25 per cent of ab initio lines in our 
1250 K list by using a Gaussian random uncertainty. This 
should represent a worst case scenario because we have re- 
moved those lines which are weaker than 1/10,000 of the 
strongest line and which are expected to exhibit the largest 
frequency /wavelength errors. We then carry out simulations 
by injecting a planetary signature into the HD 189733b spec- 
tra using our 1250 K spectrum and deconvolving firstly with 
the matching line list (case A), and the with an adjusted line 
list which models the wavelength uncertainties (case B) . The 
mismatch (i.e. case B relative to case A) leads to a planet 
which is detected with a 6.5 per cent underestimation of 
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Figure 6. Relative contrast ratio vs ad hoc deconvolution model. 
The standard model is represented by model 0. A fake planet was 
injected into the HD 189733b timeseries with 99.9 per cent confi- 
dence and recovered using our procedure. For each of models 1-4, 
the deconvolution step was carried out using the corresponding 
model. The solid line represents the recovered planetary signature 
relative to the standard model 0. The dashed lines are the 95.4, 
99 and 99.9 per cent confidence levels (bottom to top) respec- 
tively. Note how the contrast ratio is incorrectly recovered with 
models 1-4, while the relative significance increases relative to 
the standard model. 



contrast ratio and a 14.5 per cent decrease in significance 
relative to the 99.9 per cent confidence level. 



5.3 Temperature uncertainty 

The effect of using a model line list for deconvolution which 
varies from the observed spectrum in temperature alone is 
shown in Fig. 7. Here, a 1250 K planetary spectrum sig- 
nature is recovered with 99.9 per cent confidence using a 
1250 K modified BT2 line list. However both the recovered 
contrast ratio and significance change when deconvolved 
with line list temperatures which differ by ±250 K and ±500 
K from 1250 K. The effect is again relatively small for an 
underestimation of temperature (the increase in significance 
is likely due to over-weighting of strong lines and under- 
weighting of weak lines as described in §5.1) while slightly 
more significant for overestimation of temperature. In all 
instances a planetary signature is however recovered. For 
HD 189733b, the above effects alone are not sufficient to 
explain the lack of planetary signature (using our standard 
model) which is predicted at the logio(en) = -3.163 level. 
Combined wavelength and temperature mismatches should 
lead to a 99.9 per cent planetary signature appearing with 
95.4 per cent confidence at worst. A ±250K mismatch in 
model spectrum temperature results in a 20.5 per cent rela- 
tive uncertainty in the confidence of a recovered signal. 



5.4 Other sources of uncertainty 

Additional model and observational uncertainties may con- 
tribute to an incorrect estimate of the planet/star con- 
trast ratio or its relative significance. A possible source 
of error may arise from the planet ephemeris although 
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Figure 7. Significance of recovered planet as a function of tem- 
perature. A planet with 99.9% significance is simulated with a 
temperature of 1250K and recovered with 99.9% significance. By 
deconvolving the spectra with a range of temperatures, the re- 
covered contrast ratio and significance of the planet are seen to 
vary. 



this has been determined to high precision. Following our 
previo us study of the HD 18973 3 system llBarnes et al.l 
l2007bh . we adopt the ephemeris of I Winn et all (|2007l ). de- 
termined from Stromgen b and y passband observations 
(It = 2453988.80336(23) + 2.2185 733(19)). A more re- 
cent estimate of the ephemeris by lAgol et "ah (2008) us- 
ing Spitzer 8 [im observations of planetary transit yields 
T t = 2454279.436741(23) ±2. 21857503(37). Since limb dark- 
ening and starspot effects are reduced at longer wavelengths, 
this has been claimed as the most precise measurement of 
the ephemeris to date. The predicted mid-transit time for 
our observations differs by 73 sees when comparing the two 
ephemerides, a phase difference of 0.00038. The level of pre- 
cision of ephemeris observations from HD 189733b is there- 
fore now sufficient that new refinements have no measurable 
effect on the contrast or velocity amplitude of a planetary 
signal. 

A more important consideration arises from global 
re-circulation patterns in the atmosphere of HD 189733b. 
For a tidally locked planet with a static atmosphere, one 
might expect the maximum planet/star contrast ratio to 
occur at orbital phase (f> — 0.5 due to the highest ef- 
fective irradiat i on of the planet at the sub stellar point. 
iKnutson et al.l (|2007h found a difference in day and night 
side brightness temperatures of 238 K from 8 urn pho- 
tometric Spitzer light curve variations. The 1212 ±11 
K dayside temperature was found to be displaced from 
the substellar point by 16 ± 6°. Th is finding is in ac- 
cordance with 3D circulation models Showman fe Guillot] 



J2002I); ICooper fc Showman! (|2005l V, iFortnev et all (12006 
IShowman et all (|2008h . More recentlv. lKnutson et all (|2009h 
have re-analysed their 8 [im lightcurves and published 24 u.m 
lightcurves of HD 189733b. Maximum brightness is found to 
occur at phase 0.396 ± 0.022 corresponding to a shift east- 
ward of 20° - 30°of the hottest region relative to the substel- 
lar point. We have carried out a simulation to estimate the 
effect of such a shift which is not accounted for in the preced- 
ing analysis. We created an artificial planetary signal which 
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peaked 30° before secondary eclipse and recovered with a 
phase function which peaked at the same shift and also at 
secondary eclipse. The recovered planetary signature which 
did not account for the 30° shift was found to overestimate 
the contrast ratio by ~ 10 per cent. The relative change in 
significance increases by 4.7 per cent since the contrast ratio 
must be increased to optimise the \ 2 fit to the mis-aligned 
phase function. We note that this effect will be dependent 
on observational phase coverage and S/N ratio from night 
to night (i.e. shifting the phase function peak to a region of 
fewer observations or lower S/N ratio will reduce sensitiv- 

ity). 

We have assumed an effective vsini for HD 189733b of 
2.53 km s _1 which corresponds to a tidally locked planet. 
However there may be ad ditional broadening as a result of 
the re-distribution of heat. IShowman et al.l (2008) find that 
up to 3-4 km s" 1 wind speeds are responsible for the ad- 
vection of heat away from the substellar point. This shift is 
somewhat less than our resolution element of 11.99 km s -1 . 
Although the wind speeds are effectively translational (an 
east-west flow) at the 100-1000 mbar levels from which the 
2.2 [im spectrum is exp e cted to predominantly arise (see 
Fig. 4 of |Sh owman et all (2008)), we have simulated an ad- 
ditional 4 km s -1 broadening of the spectral lines. Combined 
in quadrature with the rotational broadening, we simulate 
a planetary atmosphere which possesses lines broadened by 
4.73 km s -1 rather than 2.53 km s _1 from rotational broad- 
ening alone as in the preceding sections. As expected this 
effect is also minor at a resolution of 25,000 with a 9.6 per- 
cent drop in sensitivity. 

5.5 A Semi-empirical approach - an L dwarf 
spectrum 

In addition to model uncertainties, we have carried out a 
semi-empirical examination of our ability to recover the 
spectrum of a brown dwarf which closely matches the plane- 
tary temperatu re of 1250 K. K band obs e rvations of an L3.5 - 
L4.5 spectrum (|Kirkpatrick et al.ll2000l ; iKnapp et al.ll2004l ) 
were secured by with NIRSPEC at a spectral resolution of 
22,000, covering redder wavelengths in each order than the 
HD 189733b observations. We were thus unable to use the 
L spectrum as a template which could be injected into our 
timeseries to mimic the signature of a fake "planet" . Instead, 
using the Taylor expansion scaling technique described in 
§2.3, we scaled the standard HD 189733b model spectrum to 
give the closest possible match to our observed L spectrum. 
Being an L dwarf, our spectr um exhibits significant rota- 
tion, with v sin i = 32 kms -1 (|Zapatero Osorio et a.1.1 [20061 ; 
iReiners fc Basril |2008h . The same broadening was applied 
to our standard model prior to scaling it to the L spectrum. 
Deconvolution was then carried out using: (a) the standard 
model line list on the scaled standard model spectrum and, 
(b) the standard model line list on the L spectrum. Closely 
matched deconvolved profiles are recovered in both instances 
but with a smaller equivalent width for case b. Since usini 
is matched, the resulting profiles essentially differ in their 
depths only, with the case a profile being 55 per cent deeper 
than case b profile. It is difficult to assess wavelength mis- 
match effects given the broad nature of the profile; however, 
we can attribute the 55 per cent decrease in profile strength 
to line strength mismatches. Although there may be dif- 



ferences between a L spectrum and a planetary spectrum, 
this semi-empirical approach may be taken to represent an 
upper limit to our line depth sensitivity. The line depth, 
wavelength and temperature uncertainties in §5.1, 5.2 & 5.3 
yield a 28 per cent reduction in sensitivity when combined 
in quadrature. The semi-empirical analysis result may be 
equated with this combination of effects, and is almost twice 
the modelling estimate. 



6 SUMMARY & DISCUSSION 

We have carried out a high resolution search for the 
signature of the close orbiting extrasolar giant planet, 
HD 189733b. Our signal enhancement technique enables us 
to achieve the sensitivities required to detect the dayside 
spectrum of the planet that has already been observed at a 
mean contrast ratio of F p /F* = 1 /1930 by S09 in the K band 
region of our observations. Inclusion of augmented CO2 
abundance is however not sufficient to detect the planet with 
a 99.9 percent confidence level of F p /F„ = 1/1920 (i.e. al- 
most identical to the S09 result). A tentative candidate plan- 
etary signature is found at 15 kms^ 1 greater than the ex- 
pected velocity amplitude of the planet at F p /F* — 1/4270. 
In light of the model uncertainties that have been investi- 
gated, finding a planetary signature with modified contrast 
ratio and velocity amplitude is reasonable. This prompted 
us to perform simulations in which planetary signals were 
injected at contrast ratio levels equivalent to those induced 
by contiguous absorption residuals. While these planetary 
signatures could be recovered, we found that the velocity 
amplitude may be uncertain by ±20 kms -1 , further reflect- 
ing the difficulty of reliably extracting a real signal at the 
95.4 per cent level. We note however that a planetary signa- 
ture with 99.9 per c ent confidence shoul d easily be detected, 
as demonstrated in I Barnes et al.l l|2008l ). 

Splitting the data into two wavelength regions revealed 
that the 2.21-2.36 u.m region (containing mainly H2O and 
CO opacities) yielded a candidate planetary signature with 
higher confidence. Analysing these subsetted data on a night 
by night basis however revealed that the signature was not 
stable in velocity amplitude or contrast ratio suggesting that 
it could result from a chance alignment of a number of sys- 
tematic contiguous absorption residual features at the phase 
dependent velocity position of the planet. By removing these 
features we found that the signature, close to the known 
K p = 152.6 kms -1 , disappeared. The remaining 66 per cent 
of the data did not possess the ability to recover a plan- 
etary signature at the level determined by the results of 
ISwain et al.l (|2009l ) with between 95.4 and 99 per cent con- 
fidence. Since the remaining data contained contiguous re- 
gions with levels above the mean, this may not be surprising 
as we only search for absorption signatures. It is important 
to emphasise that the tests we have carried out do not rule 
out the possibility that a true planetary signal is contained 
within the spectra. The detected candidate features may be 
partially influenced by a true planetary signature, but at the 
95.4 per cent levels, no confident claim for a detection can 
be made. 

The effects of model opacity strength uncertainties, 
wavelength uncertainties, temperature mismatch, phase 
function mismatch and velocity field/broadening uncertain- 
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ties contribute sensitivity uncertainties of 12.5, 14.5, 20.5, 
4.7 and 9.6 per cent respectively. Combining these effects 
in quadrature yields a total uncertainty in the significance 
of the result of 30 per cent. Further, if we take the semi- 
empirical 55 per cent uncertainty as an upper limit to 
our line depth, wavelength and temperature mismatches, 
the corresponding uncertainty is 56 per cent. Assuming 
that the L spectrum can provide a close match to that 
of HD 189733b, the semi-empirical result already shows 
that the model uncertainties may be significantly underes- 
timated. Hence the 99.9 percent confidence with which we 
reject a signal at the know K v could in fact be modified to 
a level with reduced significance, taking a candidate signal 
to contrast ratios that are plagued by systematic features. 

While the current generation of models can adequately 
fit broadband photometric and low resolution spectroscopic 
observations, it is clear that moving to higher resolution 
requires further model refinement. With the uncertainties 
investigated above, we can not rule out the presence of the 
planet using our technique, especially if further model uncer- 
tainties remain unaccounted for. Only further observations 
which would bring about an increase in sensitivity, or more 
precise model atmospheres could increase our chances of de- 
tecting HD 189733b. 
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